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Abstract 

We investigate the non-perturbative renormalization group behavior of the 
gauge coupling constant using a truncated form of the functional flow equation 
for the effective average action of the Yang-Mills-gravity system. We find a non- 
zero quantum gravity correction to the standard Yang-Mills beta function which 
has the same sign as the gauge boson contribution. Our results fit into the picture 
according to which Quantum Einstein Gravity (QEG) is asymptotically safe, with 
a vanishing gauge coupling constant at the non-trivial fixed point. 
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1 Introduction 



Recently a lot of efforts went into the computation of gravitational corrections to the 
beta function of the running Yang-Mills coupling constant. Robinson and Wilczek [1,2], 
in an effective field theory setting, obtained a non-zero correction at the one-loop level. It 
has the same negative sign as the familiar term already present in absence of gravity, and 
so it would render even pure abelian theories asymptotically free. After Pietrykowski [3] 
had realized that this result is gauge fixing dependent, Toms [4] reanalyzed the problem 
using the Vilkovisky-DeWitt method. In a manifestly gauge invariant as well as gauge 
fixing independent formulation of the effective action he finds that the quantum gravity 
contributions to the running charge vanish. 

All of these computations employ the dimensional regularization scheme. In [5], 
Ebert, Plefka and Rodigast pointed out that its use might be problematic since it is in- 
sensitive to quadratic divergences, and it is precisely such quadratic divergences that are 
responsible for the non-zero result obtained in [1,2]. Using a cutoff regularization instead 
they found that all gravitational quadratic divergences cancel so that there is again no 
correction to the beta function. Thereafter Tang and Wu [6] argued that the use of a 
cutoff regularization is not permissible here since it does not respect gauge invariance. 
Performing a calculation in a scheme which both retains quadratic divergences and pre- 
serves gauge invariance ("loop regularization") they obtained a non-zero gravitational 
correction to the one-loop beta function. Furthermore, Toms [7] demonstrated that, with 
a cosmological constant included, also dimensional regularization yields a non-vanishing 
gravitational correction, albeit of a different type. 

In the present paper we shall analyze the running of the gauge coupling constant in 
the framework of the Asymptotic Safety approach to quantum gravity [8-41]. Contrary 
to the calculations mentioned above it considers metric gravity not merely an effective 
but rather a fundamental quantum field theory, with the continuum limit taken at a non- 
trivial renormalization group (RG) fixed point [8-10]. Instead of perturbation theory, 
the main tool will be the gravitational average action [11] and a suitably truncated form 
of the associated functional RG equation (FRGE). Originally developed for matter field 
theories [42-47] the effective average action turned out an ideal tool for investigating the 
RG flow of Quantum Einstein Gravity (QEG) [11-38] and exploring its potential physics 
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implications [50-61]. 

The purpose of this paper is twofold: First, we develop a general framework for both 
exact and approximate ( "truncated" ) investigations of the Yang-Mills-gravity system by 
means of a gauge invariant running effective action. In particular we shall see that because 
of the semi-direct product structure of the pertinent gauge group there arises a subtlety 
as for the appropriate construction of the ghost action. Second, we use the resulting 
framework in order to find the RG flow in a simple truncation of the space of actions 
which, however, is general enough to allow for an approximate determination of the beta 
function of the Yang-Mills coupling constant in presence of quantized gravity. 

The remaining sections of this paper are organized as follows. In Section 2 we 
describe the general setup of the effective average action for a Yang-Mills field coupled to 
gravity, the corresponding FRGE in particular. Then, in Section 3, we discuss the problem 
of background covariant ghost actions and explain its solution. The explicit computation 
of the running gauge coupling is performed in Section 4, and Section 5 contains a summary 
as well as a brief discussion of various immediate applications of our results. 



Schematically the dynamics of the Yang-Mills-gravity system is governed by the 
path integral 



Here 7^ and A® are the quantum metric and the quantum gauge field, respectively, 
and S denotes the bare action. As usual, both of these fields are supposed to transform 
tensorially with respect to diffeomorphisms, 5t>- In addition, A" defines a connection 
with respect to Yang- Mills gauge transformations, 5 Y m- Denoting the vector field that 
generates the diffeomorphism by and the parameter of the Yang-Mills transformation 
by A a , we have (C v denotes the Lie derivative along v^): 



2 The Functional RG Setup 




(2.1) 




(2.2) 
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5 YM (X)Al = -d,\ a + r bc \ b Al 



(2.4) 
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From now on we will assume the Yang-Mills gauge group to be SU(N), so a runs from 
1 to iV 2 — 1, f abc are the associated structure constants. We demand S to be invariant 
under both <5 D and <5ym- 

Employing the background formalism, the dynamical fields are decomposed accord- 
ing to 

7^ = ^ + V and Al = Al + al (2.5) 

with fixed, but arbitrary background configurations g^ v and A® and fluctuations and 
[62] . Assuming a translational invariant measure, the fluctuation fields will replace the 
full quantum fields as the variables of integration in (2.1). 

There are now two possibilities to realize the gauge transformations 5d and <5ym at 
the level of the background decomposition: 

• The background gauge transformations 5® and 5yM are defined such that under 
diffeomorphisms all the fields transform tensorially, i. e. 

8*(v)$ = C V <S>, Qeig^h^A^aft (2.6) 

With respect to Yang-Mills transformations, the background gauge field transforms 
as a connection and the fluctuation transforms homogeneously: 

^m(A)4 = ~d,X a + / abc A 6 4, 5%(X)a; = f^\ b al (2.7) 

• On the other hand, we can define true gauge transformations 5% and <5yM by requir- 
ing that these shall only affect the fluctuations, not the background fields: 

0, 5g(u)V = A,fe* + V) ( 2 - 8 ) 
0, 6%(v)a; = £ v (A; + a;) (2.9) 

0, <C(AK = ~d,\ a + / abc A 6 (4 + a«) (2.10) 

The crucial idea is to choose a gauge fixing term S gi that breaks only the true gauge 
invariance but retains background gauge invariance. It gives rise to an associated ghost 
action S gh in the usual way. Introducing ghost fields and C M for the diffeomorphisms 



^m(A)A« = 
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and S a and S a for the Yang-Mills transformations, respectively, we arrive at the following 
path integral which depends parametrically on the background fields: 

Z = J Vh^Va a ^VC^VC^Vi: a Vi: a e- s[s+h > A+a] - ssf - s ^ (2.11) 

In order to set up a functional RG equation, we follow the well-known construction 
of the effective average action [42,43,47] and add a higher derivative IR cutoff term A k S 
that is quadratic in the fluctuations: 



"per 



A k s = d d x^ 9 h, v nr[gr pa h P 

+V2J d d sv^P,£)ttfM](£) (2-12) 

with k = (327rG)~5 and G denoting Newton's constant (see below). Furthermore, adding 
appropriate source terms enables us to easily compute expectation values of the quantum 
fields. With the notation 

= (Al), S« = (a-> (2.13) 
9nu = (7a«/), V = (V) ( 2 - 14 ) 

e = (n, i» = (c») (2.15) 

T a = (S a ), T a = (S a ) (2.16) 

the background decomposition (2.5) now reads 

9pu = 9^ + K u and Al = Al + al (2.17) 

The Legendre transformation of the now /c-dependent functional InZk with respect to the 
sources leads to the effective average action then [42,43,47]: 

r fe [v, a-, e, T a , T a ; ^, A;\ = r k [ g ^ ^, a«, a;, e, T a , ^1 (2. is) 

Its scale dependence is governed by the functional renormalization group equation (FRGE) 



d t T k = -STr 



(ri 2) + ^ fc (A)) _1 (^ fc (A))] (2.19) 
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where rj^ denotes the Hessian of T k with respect to all fluctuation and ghost expectation 
values, A is some suitably chosen generalized Laplacian, and t = In k is the "RG time" . 

Since we deal with two gauge invariances, a remark concerning the notation is in 
order: We write D = d + T and V = d + A for the covariant derivatives that are 
constructed by means of = \g pa {d v g atl + d^g av - d a g^ u ) and respectively; the 
covariant derivative containing both of these connections is denoted by V = d + A + T. 
By adding a bar, we denote their analogues evaluated on the background configurations. 

3 On the Construction of 

Background Gauge Invariant Ghost Actions 

3.1 Motivation 

Since we have to fix two gauge invariances by two gauge conditions, F^(h^ u ] g^ U) A°) 
and G a (a^; g^, A®) for the diffeomorphisms and the SU(N) transformations, respectively, 
the associated Faddeev-Popov operator will in general consist of four components. (Here 
we have already assumed that the diffeomorphism and the SU (N) gauge condition only 
involve the metric and the gauge field fluctuation separately.) The corresponding classical 
ghost action will then be of the form 

S gh [h, a, C,C, S, Z;g,A] = -J d d x yfi (V 1 C^^-S$(C) V+ 

^ (3-1) 
+ k" 1 C^^^m(S) V + g S a ^^(C)< + g S^ym^K) 

Since we are going to neglect renormalization effects in the ghost sector, the evolution 
equation for T k will contain only the classical ghost action but with the quantum ghost 
fields replaced by their vacuum expectation values, and the full classical fields g^ u and A® 
identified with their background configurations g^ u and A a ^ respectively; stated differently, 
in this class of approximations the fluctuations and a® can be set to zero in even 
before is computed. 

Looking at the T-£ part of the resulting ghost action we encounter a serious problem: 
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Employing G a (a; g,A) = V^a a as gauge condition, it is given by 



(S gh [£, T; g, A]) n = -J d d x ^ [T a V^ b (i"d p A b v + {d^)A b p )] (3.2) 



Here the covariant background derivative T>^ acts on the ordinary Lie derivative of 
an SU(N) background connection; therefore, this part of the ghost action is not 5yivr 
invariant. In order to resolve this issue, it is useful to look at it from a more abstract 
point of view. 

3.2 Ward operators and their algebra 

To begin with, we consider an arbitrary functional F^^, A^, C M , C M , S a , S a ] of the 
dynamical fields at hand. At this stage the splitting into fluctuations and background 
configurations has not yet been performed. An infinitesimal gauge transformation of F, 
considered a scalar functional of its arguments, consists of a diffeomorphism along v M and 
an SU(N) transformation with parameters A a . It can be implemented as 



^[7 + M^7 + <WA)t, A + S D (v)A + 5 YM (\)A, C + 5 D (v)C + 5 YM (\)C, 
C + 5 D (v)C + (5 YM (A)C, E + 5 D (v)Z + 5 Y m(A)S, S + 5 D ( V )S + <5 YM (A)£] 

(3.3) 

= F[ 7 , A, C, C, S, S] - V^(u)F[ 7 , A, C, C, S, S] - W YM (A)F[ 7 , A, C, C, S, S] 
+ C(v 2 ,t;A,A 2 ) 



with the corresponding Ward operators generating diffeomorphisms, 



W D (v) = - f d d x (s D (v)^ u (x)-^- 



+ 6 D (v)Al(x) 
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+ 5 D (v)C»(x) 



5 



+ 5 D (v)C l ,{x) 



5 



(3.4) 



5C»(x) 



SC^x) 




6 



and Yang-Mills gauge transformations: 

d d x (£ YM (A) 7 ^(x) — ^ + 5 YM (A)^(x)^^y 

+ ^M(A)^(,)^ + 5YM (A)C,(,)^- y (3.5) 

+ 



*vm(A) S »(,)^^ ym (A)S"(x)_) 



In these integrals the measure factor ^/det(7 M!/ ) cancels against a similar one which would 
render the functional derivatives tensorial. We define the set of all invariant functionals 
by 

F im = {F | W D (v)F = A W YM (A)F = 0V«", A a } (3.6) 

Computing the algebra of these operators leads to 

[Wd^O.WdW] = W D (K,^]) (3.7) 
[Wym(Ai),W Y m(A 2 )] = W Y m(/A!A 2 ) (3.8) 
[m(v),W YM (\)\ = Wym(C v X) (3.9) 

where [^1,^2] denotes the Lie bracket and {f\i\2) a = f abc X\^2 ■ This algebra implies that 
the total group of gauge transformations, G, has the structure of a semi-direct product 
of the spacetime diffeomorphisms Diff and the local Yang-Mills transformations SU(N) loc , 
with the latter forming the invariant subalgebra: G = Diff x SU(N) loc . Whereas the first 
two relations (3.7), (3.8) represent the well-known composition laws of diffeomorphisms 
and Yang-Mills gauge transformations, the third relation (3.9) lies at the heart of our 
problem: diffeomorphisms and local gauge transformations do not commute. Instead, 
they close on the Lie derivative of the gauge parameter. In particular, this implies that 
diffeomorphisms do not map SU (N) tensors onto SU (N) tensors. 

3.3 Modified diffeomorphisms 

What is called for is an SU(N) covariantization of the ordinary Lie derivative. This 
is tantamount to a different parametrization of G that makes the mixed commutator 
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vanish. This can be achieved by defining new diffeomorphisms which include a SU(N) loc 
transformation with parameter A a = A^v M [63] : 

Wv(v) = Wd(v) + Wym(A ■ v) (3.10) 

Loosely speaking, this amounts to shifting a certain v -dependent part of SU(N) loc into the 
diffeomorphism sector. The notion of an "invariant functional" remains unchanged: 

•Finv = {F | Wd(v)F = A W YM (A)F = 0V/, X a } (3.11) 

The algebra relations receive extra contributions now since the Ward operators act on the 
field dependent parameters of the transformations as well. This leads to an algebra of the 
desired form (with (viv% ■ F) a = v^v^F^): 

[VVdK),VVd(^)] = m>([vi,v 2 ])-W YM (v 1 v 2 -F) (3.12) 

[Wym(Ai),W ym (A 2 )] = W YM (/A 1 A 2 ) (3.13) 

[m>(v),W Y M(\)\ = (3.14) 

If we now split the dynamical fields into fluctuations and background configurations, 
we have to decide whether the field dependent transformation parameter in (3.10) should 
contain the full or the background gauge field only. Since, as already mentioned, the 
metric and gauge field fluctuations do not enter the final form of anyhow, we may 
safely opt for the latter already at this point: 

M> G {v) = W^ G (v) + W^(A ■ v) (3.15) 

According to the usual distinction between true and background gauge transforma- 
tions, we now have to consider two classes of Ward operators as well. In the background 
case, the algebraic relations simply carry over, so we have for the Ward operators gener- 



ating background gauge transformations: 

[yVE(vi)M(v2)} = yVE([vi, v 2 \) - >Vym(^2 • F) (3.16) 

[W?m(Ai),W? m (A 2 )] = Wft^Aa) (3.17) 

[WE^MuW] = (3.18) 
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As for the true gauge transformations, the merely background field dependent pa- 
rameter A a = A^v^ of the compensating SU(N) loc transformation is not subject to true 
gauge transformations. The algebra that generates these transformations can be easily 
computed therefore by taking advantage of the linearity of commutators. We obtain 



[W^m(Ai),W« m (A 2 )] 
[WgW,W£ M (A)] 



yvE([vi,v 2 ]) + yv^ M (v 1 v 2 -F) 

W£*(/AiA 2 ) 
W^ M (v • VA) 



(3.19) 
(3.20) 
(3.21) 



By writing W we distinguish these Ward operators from the modified one that was defined 
with respect to the undecomposed fields; in addition, their notation shall remind us of 
the fact that the gauge field entered their definition only via its background component 
A. 

Thus the actual theory space on which we can define an RG flow consists of the 
functionals F[h, a, f , £, T, T; g, A] = F[g, g, A, A, £, £, T, T] in 



Kv = {F I W»F = A W* M (\)F = 0V/, A a } 



(3.22) 



Finally, we return to our starting point and compute the Y-£ part of the ghost action, 
now with the original true diffeomorphism 5§ replaced by its modified counterpart, 5*g'- 

{S gh [^f;g,A]) n = - J d d xV§[r a V^ b (5§(Oa b u )]\ a=0 



-J d d x ^[r a V» ab (ed p Al + (d u e)A 

- (d u (A b p e) + f bcd A b A d p e 
-J d d x^[-r a v» ab F b p e 



(3.23) 



This action is obviously invariant under SU(N) loc and background diffeomorphisms for 
background tensorial ghost expectation values. 
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4 The Running Gauge Coupling 



In this section we explicitly evaluate the FRGE on a truncated theory space which 
is general enough to allow for an approximate determination of the beta function for the 
scale dependent Yang-Mills coupling <?ym(&)- Our truncation is given by the following 
ansatz: 

F k [g, g, A, A, £, I T, T] = Tf H [g} + r Y k M [g, A] + Tf[g -g,A-A;g, A] 

(4.1) 

+S gh [g-g,A-A,dT,t;g,A] 

Here 

r™[g] = 2n 2 Z N {k) J d d x y/g (-R(g) + 2X(k)) (4.2) 

is a /c-dependent form of the Einstein-Hilbert action. The corresponding dimensionful 
running parameters are the cosmological constant \{k) and Newton's constant G(k) = 
G/ZN(k) where G is a fixed reference value. Furthermore, 

rn<?, A] = §^ / d d x Vgg^F^ (4.3) 

is the standard second-order Yang-Mills action, with a A;-dependent prefactor Zp(k) 
though. Hence the (dimensionful, except in d=4) running gauge coupling is ^ym(^) — 
3ym^f(^)~^ 2 with some constant gxu- Finally, 

rf [g -g,A-A;g,A\= f d d x ^ (^9^K + f^G'G") (4.4) 

implements the gauge fixing conditions for the diffeomorphisms, F M , and the SU (N) gauge 
transformations, G a . As in [11] and [43] we factored out the wave function renormaliza- 
tions Z N and Z F from the gauge fixing parameters and o:ym, respectively. In principle 
the latter are still /c-dependent but we shall neglect their running here. In fact, later on we 
set o:d = «ym = 1- Our choice for the gauge conditions complies with the requirements 
discussed in the previous section: 

F,(h; g) = V2k (sfa^D^ - \g aP D^j h af3 (4.5) 
G a (a;g,A) = g~Y u V,K (4.6) 



10 



The resulting ghost action reads, with a° 7^ and h pv ^ still, 



Sgh[h,a,£,£,r,r;g,A] = 

- J d d x J=g (v^ {r P 9 aX Dx (g P vD a + 9av D p ) - g p °r X D x g au D p ) C 

+ t a g^V p (F; u e + ?d p al + {d^)a a p + f abc A b p ^a c u ) + t a {g^5 ab V p V p ) T 6 ) (4.7) 

It can be checked that S g h of eq. (4.7) is invariant under background gauge transfor- 
mations: WyM'S'gh = = WpiSgh. While this is true even for non- vanishing fluctuations 
h and a, in the present calculation we shall need S s h only for h = = a. 

At this point a remark concerning the expected reliability of this truncation ansatz 
might be in order. As for its gravitational part, all generalizations of the Einstein-Hilbert 
truncation explored during the past decade did not change the qualitative picture it gives 
rise to, at least close to the non-Gaussian fixed point. In the Yang-Mills sector we re- 
tained only the first monomial of a systematic derivative expansion. From the analysis 
in [43] without gravity we know that this truncation is not only sufficient to reproduce 
one-loop perturbation theory exactly, but even approximates the two-loop result for the 
beta-function with a small error of a few percent. Therefore we may expect that this trun- 
cation, too, is perfectly sufficient as long as k is sufficiently large (well above dynamically 
generated "confinement" scales, say). 

When we insert the truncation ansatz (4.1) into the exact FRGE (2.19) the super- 
trace decomposes into a "bosonic" and a ghost contribution: 



d f T 



U k 



Tr 



Tr 



o(2) , -pgh/o—l o(2) \ 
^gh + X-k l^gh °gh ) 



(4.8) 



Here f k = Tf H + Tj M + T| f is the bosonic part of the action and f {2) is its Hessian. The 
coarse graining operators in (4.8) have the structure 



U k (x) = Z k k 2 R^\x/k 2 ) nf{x) = Zfk 2 R^\x/k 2 ) 



(4.9) 



where R^\y) is a "shape function" continuously interpolating between R^(0) = 1 and 
lim R(°\y) = 0. The constants Z k and Zf 1 are matrices in field space. They will be 

y— >oo 

adjusted in such a way that if in a certain mode has the inverse propagator ( k p 2 it 
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becomes ( k (p 2 + k 2 R^) when we add 1Z k to rj^. As we shall see, this requirement is 



met if Z k and Zf h have the following block structure in (h, a, £, £, Y, T)-space: 



( Z k)- h h 

{z?) 



= 2 



(Zk)- 



afib 



Z F {k) 



9 2 



5 ab 5» 



ab 



(4.10) 



Note that Zf 1 is actually /c-independent. 

In setting up eq. (4.8) we followed ref. [43] and opted for the complete Hessian 



operator to play the role of A. 1 More precisely, we set A 



Z- x Yf and A 



7 -lo(2) 



in the (h, a)- and the ghost-sectors, respectively. The multiplication by the inverse Z 
matrices brings A closer to an ordinary (covariant) Laplacian; symbolically, if T k = 
—(kd 2 + • • • , Z k = ( k , we employ A = — d 2 + ■ ■ ■ rather than A = —( k d 2 + • • • . 

The most complicated ingredient needed in order to evaluate the traces in the FRGE 

° (2) 

is the Hessian of the bosonic action Fy, i. e. the matrix of its second functional derivatives 
with respect to the dynamical fields (h,a), or equivalently (g,A), at fixed backgrounds 
(g, A). This Hessian is most transparently displayed by means of the associated quadratic 
form ad which appears in the expansion 

f k [g + h, A + a, g, A] = f k [g, A, g, A] + 0(h, a) + r^^, a; g, A] + 0({h, a} 3 ) (4.11) 



Explicitly, r^ uad is the sum of the following terms which reflect the block structure of 



(2) 



1 Recently this kind of cutoff has been referred to as "spectrally adjusted" [46,48,64] or as of "type 
III" [23]. 
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in (h, a)-space: 

+ (l - ^ v ,D p D„ + vC \f^ h* (4.12) 
( rr") = -M- U d x^a a J- 5 ab 5^V p V, + 2g^f abc F c pv + 8 ab gp>R m + 

V / aa 1 g yM J \ 

+ ( 1 - — j S^gtPVpVr, j a 6r? (4.13) 

(4.14) 

(rr d ) K = (r quad )„ (4.15) 

\ /ah 

The above quadratic functionals contain the kernels 

KX 'vC = \ ( 5 v 5 C + 5 v 5 C - 9 X %<) (4-16) 

f** = \ (& X A + <^C - 9 X %<) {R ~ 2A) + gtRrt - 5^ c - (4.17) 

(j9 x *9 p '9nC - " \9 XP f a 9vt + <«^) (4-18) 



+ 2 {Sffig^g* - Sffig*?' + 26ffig" x 9 m ') (4-19) 

Using these formulae it can be checked that the Z-factors (4.10) are correctly chosen. 
Since we are not going to extract any "extra" background field dependence [37] we may 
set g^ = g^ and A a p = A® after having found the Hessian. 

The truncation contains three running couplings, <?ym(&), G{k) and X(k). Their 
beta functions can be found from the FRGE (4.8) by "projecting out" the corresponding 
invariants in the derivative expansion of the traces and equating them to the corresponding 
field monomials on the LHS of the flow equation. The resulting system of differential 
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equations becomes autonomous if we employ the dimensionless counterparts of <?ym, G 
and A respectively: 

L.d-2 

g 2 YM (k) = k d - 4 Z F \k)g 2 YM , 9(k)^^^-^ = k d - 2 G(k), X(k) = k~ 2 \(k) (4.20) 

In terms of these variables the three coupled RG equations have the structure 

d tg Y M = Pym = (d-A + r] F ) g YM 

d t g = P g = (d-2 + m )g (4.21) 
d t X = Px 

Here we introduced the anomalous dimensions related to the Yang-Mills and the gravita- 
tional field, respectively: 

Vf = -d t In Z F , VN = -d t In Z N (4.22) 

In this paper we are only interested in the gravitationally corrected Yang-Mills 
beta function f3 YM - Therefore it is sufficient to extract the F^-term from the derivative 
expansion of the traces. For identifying this monomial and reading off its prefactor we 
may insert any metric. We shall employ the most convenient choice, g^ = g^ = 5^. 
Furthermore, we set ctD = «ym = 1 from now on. The remaining calculation is in principle 
straightforward, but rather lengthy. One has to expand the traces up to terms with two 
fields A a (x) and two derivatives acting on them. Because of the built-in background gauge 
invariance those terms should combine to F® u F aflu . As a check we verified that this indeed 
happens. 

Let us now discuss the result. Here we specialize for d = 4 spacetime dimensions; for 
general d the reader is referred to the Appendix. We present three different formulae for 
i]f', they differ with respect to the degree of "RG improvement" they take into account. 

To start with, we "switch off" all RG improvements. This means that we discard all 
terms in d t TZk on the RHS of the flow equation where d t hits either a ^-factor or the rj^ 
in the argument of TZk- In this way the evaluation of the FRGE amounts to a one- loop 
calculation, with a non-standard regulator though. We find 

VF = -^g<S>[(0)-^LNg YM (4.23) 
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so that 



%ym = 99YM®\(fy 
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Ng$ M 



(4.24) 



24tt 2 



Here $J(0) is one of the integrals which were encountered in the pure gravity calculation 
[11] already: 



The second contribution on the RHS of (4.24) is the familiar "asymptotic freedom" term 
due to the self-interaction of the gauge bosons, while the first one, due to the virtual 
gravitons, is new. 

Several comments are in order here. 

(1) The gravitational correction is manifestly cutoff scheme dependent, i.e. it depends, 
via $i(0), on the shape function R(°\ However, for any admissable choice of the 
constant &\(0) is positive. As a result, the gravity term has a qualitatively similar impact 
on <?ym(&) as the gauge boson loops, namely to drive ^ym(^) smaller at larger k. It tends 
to speed up the approach of asymptotic freedom. 

For the exponential cutoff R^(y) = y/(e v — 1), for instance, one finds $}(0) = vr 2 /6, 
while the "optimized" one [49], R (0) (y) = (1 - y)9(l - y), yields $}(0) = 1. 

(2) The gravitational correction, in perturbative language, originates from a quadratic 
divergence or, in FRGE language, a quadratic running with k. For this reason its scheme 
dependence is by no means surprising or alarming. Rather, it is the usual situation which 
is always encountered when the effective average action is applied to matter theories 
with a quadratic running of parameters, masses, say. However, one should note that 
the couplings in Tj, as such are not observable or "physical" quantities. The latter must 
be R^ -independent. This independence comes about by a compensation of the scheme 
dependence among different running couplings. (In truncations this compensation might 
not be perfect.) In general there will be compensations between effective propagators and 
vertices, for instance. Analogous remarks apply to the gauge fixing dependence. 

(3) The beta function for g\ M depends on all three couplings, (jy M , g and A. In the approx- 
imation of (4.24) it happens to be independent of A, but it does depend on g{k) = k 2 G(k), 
the dimensionless Newton constant. Hence the differential equation for g YM cannot be 




(n > 0) 



(4.25) 
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solved in isolation. In principle the full system (4.21) should be considered, and this 
would include the backreaction of the matter fields on the running of the gravitational 
parameters g and A. We shall not study this backreaction in the present paper. Instead, 
let us assume that the complete RG trajectory k \— > (gyyiik), g(k), X(k)) admits a clas- 
sical regime [58] in which Newton's constant does not run appreciably so that we may 
approximate 

G(k) w G = const, g(k) = G k 2 (4.26) 

This approximation, implicitly, has been made in all perturbative studies [1-7]. With 
(4.26), for an abelian field (N = 0), say, 

d t g YM = --$ 1 1 (0)G k 2 g YM (4.27) 

71 

Incidentally this beta function has the same general structure as the result by Robinson 
and Wilczek [1]; it is proportional to Gog YM and depends explicitly on the energy scale k. 
Its /c 2 -dependence indicates that the underlying quantum effect is related to a quadratic 
divergence. 

Eq. (4.27) is easily solved: g YM (k) = g YM (0) ■ exp (— cuym^/^pi) 2 )- Here wym = 
3^\(0) /it and m P1 = G Q 1 ^ 2 is the (ordinary, constant) Planck mass. To first order in the 
fc/mpi-expansion we get 

9v M (fc) = <?ym(0) [1 - Lo YM (k/m Pl ) 2 + 0(k 4 /m Pi }} (4.28) 

We note that to leading order Newton's constant itself [11] has an analogous scale depen- 
dence, including the sign of the correction: G{k) = Go [1 — u(k/m P \) 2 + •••]. 
(4) In order to illustrate how the above result fits into the asymptotic safety picture of 
Quantum Einstein Gravity [8,21,41] we consider a free Maxwell field again. It is known 
that, in the Einstein-Hilbert truncation, the RG flow of the average action possesses a 
non-Gaussian fixed point for the two gravitational couplings, (g*, A*), both in pure gravity 
[11,13,18] and in presence of a free Maxwell field [22]. At this fixed point the dimensionless 
Newton constant equals a positive constant, g{k) = g*, while the dimensionful one runs 
to zero quadratically: G(k) = g*/k 2 — > for k — > oo. In this regime, 

d t g 2 YM = -l${(0) g* g 2 YM (4.29) 
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The solution to this equation reads 

g 2 YM {k)<xk- e ™, e YM = -<f$i(0) (4.30) 

7T 

At the fixed point the gauge coupling approaches zero according to a power law with a 
critical exponent Oym 5 a positive number of order unity. 2 Thus the total system has a 
non-trivial fixed point of the form (<? Y m — 0, <?* > 0, A* > 0). Obviously the approach of 
9ym = is much faster than without gravity where <7ym(&) ck l/ln(/c). Note that qym 
is a relevant parameter, it grows when k is lowered, hence it contributes one unit to the 
dimensionality of the fixed point's UV critical manifold. 

Next we present the results for r\p with the RG improvements included. In a first 
step we retain only the terms which arise when d t hits the Zfc-factors in TZk- Those terms 
are proportional to r\p and 77 at, respectively. As now r\p appears also on the RHS of the 
RG equation we obtain an implicit equation for it. Its solution reads 

_ _ -U^M-^g^-l m \g 
l-i^i(0)-d^<?YM-fA 5 

In this approximation r\p depends not only on Newton's but also on the cosmo logical 
constant. Eq. (4.31) resums terms of arbitrary order both in gyyi and g; it generalizes 
a known result [43] for pure Yang-Mills theory. If we go on and include also the terms 
coming from the scale derivative of rj^ in the argument of TZ^ we are led to 

"f 9 - ^Ng YM + ±r, N g (3I1 - 2A - 2A 2 £j) - ±g A(2A + 
Vf = jz ^ ~ (4-32) 



The eqs. (4.31) and (4.32) contain the following integrals involving R^>: 
$p( w ) = —— / dz z n - 1 - ^— ^ (n > 0) 

m w ) = (l + w)- p (n = 0) 



m w ) = —— / dz z n - 1 - ",\ r- (n > 0) 

nK ! r(n)7 [z + R(°)(z) +w]p V ; 

& ( w ) = {U (n = 0) 
ov 1 (1 + w)p 



(4.33) 



2 One cannot easily extract the precise numerical value of 6ym from existing calculations since the 
determination of g* in the Einstein-Maxwell system in ref. [22] employs a cutoff different from the present 
one. 
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The threshold functions ^(w) appeared already in the Einstein-Hilbert truncation of 
pure gravity [11]. In eq. (4.32) we abbreviated = $£(0) and = $£(0), respectively. 
What is new in (4.32) is the occurrence of the complete beta function of the cosmological 
constant, (3\. It originates from the differentiation of the A-term contained in T k . 

5 Discussion and Conclusion 

In the recent literature on the gravitational corrections to the Yang-Mills beta func- 
tion [1-7] there has been a certain amount of confusion as some of the computations do 
get a non-zero result while others don't. However, we believe that different calculations 
have no reason to yield the same result unless they agree on virtually all details of the 
regularization and renormalization procedure. The quantum effects of interest are related 
to quadratic divergences (or a /c 2 -running) , and so we should not expect the same high 
degree of universality as in the case of the familiar gauge boson contribution which is 
related to a logarithmic divergence. 

In the present paper we computed the beta function for <?ym(&), defined as a coef- 
ficient in the derivative expansion of the effective average action. This approach has two 
features which are essential here: First, it retains all quadratic divergences (as opposed 
to dimensional regularization, say), and second, by the background field technique, the 
regularization (the cutoff IZk) preserves gauge invariance. 3 In this setting, we do get a 
non-zero gravitational correction. This correction is scheme and gauge fixing dependent 
but, as we explained, this is by no means unexpected but rather the usual situation. When 
observable quantities are computed from the scheme and gauge fixing dependences will 
cancel among the different running couplings involved. 

Among the perturbative calculations only the one by Tang and Wu [6] is directly 
comparable to ours. They employ a regulator which retains quadratic divergences and 
treats them in a gauge invariant manner. It is gratifying to see that they, too, get a 
non-zero gravitational correction which has the same structure as ours when we omit the 
RG improvements. 

As a first application of our results we mention that, by a standard argument, knowl- 

3 As in all traditional applications of the average action, r& is not gauge fixing independent, though, 
i. e. we do not use the Vilkovisky-DcWitt method. 
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edge about the /c-dependence of wave function normalization constants such as Zp{k) can 
be used in order to deduce information about the related fully dressed propagator implied 
by T = r fc=0 . In the case at hand the running inverse propagator of the gauge field, 
on a flat background, has the form Zp{k)p 2 oc g^i^p 2 . At high momenta, if there is 
no other relevant physical cutoff scale but the momentum itself, the dressed propagator 
D{jp) obtains by setting k = \p\, whence Dip)' 1 oc g Y 1ii\p\)p 2 ■ For the example of eq. 
(4.28), for instance, this leads us to expect that the photon propagator gets modified by 
a p 4 -term when p approaches the Planck scale: 

Dip)' 1 =p 2 + uj YM p 4 /m 2 Pl + Oip (i /ml >l ) (5.1) 
Likewise the fixed point running of (4.30) implies the following behavior for p 2 — > oo: 

Dip) oc i/ p 2 (n-e Y M/2) ( 5 2 ) 

As ©ym is positive the gauge field propagator falls off faster than 1/p 2 , thanks to the 
quantum gravity corrections. In fact, the same argument when applied to the graviton 
propagator leads to a l/p 4 -behavior for p 2 — > oo [13]. The asymptotic propagator (5.2) 
suggests that the quantum gravity corrections improve the finiteness properties of the 
matter field theory, and this precisely fits into the picture of asymptotic safety. It is also 
interesting to note that (5.2) leads to a modified static electromagnetic potential -A (r) of 
a classical point charge. 4 The 3-dimensional Fourier transform of (5.2) yields the potential 
A (r) oc r 0YM_1 which, if ©ym is large enough, could even be regular at r = 0. This makes 
it obvious that the gravity induced running of the gauge coupling is closely related to the 
old problem of divergent self energies. 



Treating the source dynamically, also form factor effects need to be included [65] 
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A Appendix 



In this appendix we display the equation for tjf in d spacetime dimensions. With 
all RG improvements included, it assumes the form: 

(47r) d /V = -32tt<? A X| _ 10 (-2A) + Bd- d / 2 Xi _ 1)0 (-2Ad) 
+ C(2A)- 1 (x|, (-2A) - x- >0 (0)) 

+ D {^/ 2 Xf+li2 (-2Arf) - Xi+1)2 (0) - \ (d- d / 2 X ^{-2\d) - Xi ,M) 



-9ymN 



26 -d 



X|- 2 ,o(°) 



Ax i _ 1)0 (-2A)+ J Bd- d / 2 X|_ 1;0 (-2Ad) 



+ C ((2A)- 1 x |i0 (-2A) - (2A)- 2 (x- +1>0 (-2A) - x- +1>0 (0))) 



d- d ' 2 



d -d/2+l 



+ 16tt7] Fg 



j—[Xl+i,2(- 2Xd "> + — 2 — i(-2Ad) - ^— j;Xf+i, 2 (0) 
C (-(2A)- 1 x f , (0) + (2A)~ 2 (x f+1;0 (-2A) - x< +1>0 (0))) 



d -d/2+l 



d 



-^-X| + i, 2 (-2Ad) - -X|,i(0) + ^3yX-+i )2 (0) 



,'24 -(f. 
+ VFg Y M N — n — X#_2,o(°) 



6 



- 167t(?7f - r] N )g 



A^ (-2\) + Bd- d / 2 ^ (-2\d) 



+ C(2A)- 1 (x| +1 ,o(-2A)-X| +1 ,o(0) 
+ ^(^2) {^ +11 (0) - ^ |+lil (-2Arf)^/ 2 } 



+ 32ng(2\ + p x 



(-2X) + Bd- d ' 2+ %_ lfi (-2\d) 



+ C ((2A)- 1 | |;0 (-2A) - (2A)- 2 (|| +1 ,o(-2A) - *< +1>0 (0))) 

~ 2Dd |_^ |+1>2 (-2Ad) - _^ +1>2 (0) + -cT^ (1 (-2Ad) 



(A.l) 
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Eq. (A.l) contains the following functions of d: 



(d + 2)(d 2 -9d+12) (rf-4)(rf-6) 
4d ' 2(d-2) 



^_ 4(3rf 2 -10rf + 4) 4(rf-4) 2 
^ ~~ tt - ; — ^ i/ = 



(A.2) 



d(d-2) d(d-2) 2 (d-l) 
Furthermore we introduced the following new threshold functions: 

{z - wY' 1 BS°\z) -zRW(z) 



Xn,m(w,v) = — -/ dz 





1 






>- 


to) 




1 




r< 


[n- 


to) 




1 





"OO 



' w 



(z + v) m z + R<®(z) 



00 (z — wY' 1 R(°Hz) 
Xnm(w,v) = —^ -/ dz^- J - , v n , . (n>m>0) (A.3) 



A "' mv ' ' V{n-m)J w (z + v) m z + R(°)(z) 

When used in Eq. (A.l) with one argument only, the second argument of the x functions, 
v, is set to v — — Jzf implicitly, for clarity purposes. 

For n G N and to = the x integrals are independent of v and can be expanded in 
w into a finite sum of $ integrals, 

fc=0 

and analogously for x and For d = 4 only x integrals that are expandable in this way 
occur. 

For to > a non-integrable singularity might occur in the integration interval (de- 
pending on the signs of w and v). In this case the definition of the x integrals is meant 
merely as an abbreviation, in the following sense. The integrands of the x integrals em- 
braced by curly brackets in (A.l) should be summed first and only then be integrated 
over the intersection of their integration intervals. This removes the apparent singularity 
and we are left with a sum of finite integrals. 

Eq. (A.l) can be reduced to the corresponding 1-loop result by omitting from its 
RHS all terms containing to?, r] N , or any x integral. The latter arise from the differentia- 

^ (2) 

tion of rjj. ; in the arguments of 1Z k . 
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